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Abstract 

Transverse momentum distribution of gauge bosons and is 
" ^ ■ calculated using unintegrated parton distributions obtained by solving 

CD . leading order CCFM equation(s) in the impact parameter space. For 

\ illustration we compare the results of the fixed-order collinear, soft- 

gluon resummation and CCFM approaches. A parameter of a nonper- 
turbative form factor is adjusted to the DO data. In contrast to the 
collinear approach, the leading order calculation based on the unin- 
tegrated CCFM parton distributions generates transverse momentum 
distributions of gauge bosons, which are almost consistent with ex- 
perimental data provided the parameter of the nonperturbative form 
factor is appropriately adjusted. The W~^-W~ asymmetry is calcu- 
li ■ lated as a function of the boson rapidity and transverse momentum. 

Predictions are given for RHIC and LHC. 

• 1— I - 
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<^ : 1 Introduction 

The production of gauge bosons at hadron colliders provides fundamental 
tests of perturbative QCD. The understanding of transverse momentum dis- 
tribution is of particular interest in this context. The data from the Teva- 
tron collider cover regions of both small and large gauge-boson transverse 
momenta {pt) At large pt the standard fixed-order pQCD approach 

is adequate. At small pt of vector bosons {pt ^ My) one encounters large 
logarithms of My /pi. The leading term is of the form 



^^^In^ (1) 
dpi p] p1 
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The logarithms can be resummed to give a Sudakov form factor. In 
recent years the resummation approach became the standard in calculating 
transverse momentum distribution of gauge-bosons [U IHl El 13 El IH]- The 
most often used is the formalism proposed by Collins, Soper and Sterman 
13]. In the following we shall refer to it as CSS or "the b-space resummation 
approach". The b-space resummation approach is subjected to an ambiguity 
in setting the matching condition between the resummed and fixed-order 
results. Different prescriptions have been used throughout the literature. 

The coUinear parton approach was used by a few specialized group to cal- 
culated many high-energy observables. Not all of them can be fully explained 
in the fixed-order coUinear parton model. The fct-factorization approach, us- 
ing unintegrated parton distributions, is an alternative approach. In recent 
years there were many calculations of specific observables using different sets 
of unintegrated parton distributions. Mainly observables sensitive to uninte- 
grated gluon distributions have been discussed, like the dijet or heavy-quark 
pair production in high-energy electron-proton scattering at HERA or the 
production of heavy quark pairs in hadron-hadron collisions at high energy. 
The existing unintegrated gluon distribution differ one from another. In some 
approaches (models) the unintegrated gluon distributions are limited to the 
small-x region by construction. The unintegrated (anti) quark distributions 
are very sparse and often not available. Clearly, the program of a system- 
atic study of the world data in terms of unintegrated parton distributions 
(uPDFs) is at its beginning. 

Gluon radiation along the ladder generates transverse momentum of par- 
tons through the recoil effect. Evolution of the unintegrated distributions is 
described then by the so-called CCFM equation [Tn]. The CCFM equation 
has a nice virtue that it interpolates between the BFKL evolution at very 
small values of x and the DGLAP evolution at intermediate and large values 
of X. For many practical applications it is very useful to solve the CCFM 
equation by means of a Monte Carlo method The CCFM approach turns 
out very succesful in resolving some high-energy puzzles jTHl HZl UHl UH] • In 
the so-called b-representation the CCFM equation for b-dependent uninte- 
grated distributions takes a simple diagonal form |llj. This approach can 
be generahzed to include in addition to unintegrated gluon distributions also 
unintegrated (anti)quark distributions jT3]. 

The similarities and differences betwen the standard soft gluon resum- 
mation and the CCFM resummation have been discussed recently on the 
example of the gg H fusion mechanism It has been shown how the 
soft gluon resummation formulae can be obtained as the result of an ap- 
proximate treatment of the solution of the CCFM equation in the so-called 
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b-represent at ion . 

We shall show that in contrast to the coUinear approach the use of uPDFs 
leads to a non-5-like transverse momentum distribution of the produced 
gauge bosons already in the leading order for matrix elements. It is our aim 
here to calculate the transverse momentum distribution of the charged gauge 
bosons based on the unintegrated parton distribution originating from the 
numerical solution of the CCFM equation. This will be a first phenomenolog- 
ical test of the CCFM parton distributions from ^Hj- In addition, we wish 
to concentrate on the — W~ asymmetry, especially for proton-proton 
collisions at RHIC and LHC. We shall show that this asymmetry strongly 
depends on the d — u asymmetry of parton distributions, providing a new 
interesting tool to limit the latter somewhat better. 



2 Formalism 

In the fixed-order perturbative treatment the situation is as follows. In the 0- 
th order coUinear calculation the produced gauge bosons carry no transverse 



momentum, i.e. 



'^^ WSr^Ur |2 

2^ IV 



dyd?pt 

qq' 



For brevity, above we have introduced = . The quark/ ant iquark mo- 
mentum distributions are evaluated at momentum fractions Xi.2 = exp(±?/). 

A nonzero transverse momentum of gauge bosons is obtained only in 
the next-to-leading order calculation. Then the gauge bosons with nonzero 
transverse momenta are produced in the weak-QCD annihilation or Comp- 
ton processes in association with a gluon (annihilation) or quark/antiquark 
(Compton). The cross section associated with the 2 — >• 2 subprocess (1+2 
W+4) is given by the standard parton formula 

where the matrix elements squared, averaged over spins and colors of in- 
coming partons and summed over spins and colors of final particles read: 

y tu 



qq' 
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for the annihilation and 



E W,,-.w,? = rcas{f^l)V2GFM'M,'\' 1 ^' + ^' + ^^^^ (5) 

su 



qq' 



for the weak Compton subprocess. In the equations above Vqqi are elements 
of the Cabibbo-Kobayashi-Maskawa matrix. The parton distributions are 
evaluated at xi_2 ~ ■ (mj^vK exp(±?/i) +ptexp(±y2)) where the transverse 

mass mt^w = \/M^ + p1- 

The poles at t = or -u = cause the theoretical cross section to diverge 
as — and exceeds the experimental data at < 5 GeV (see Fig^. 
Virtual corrections to qq' W contribute only at pi = [da^ /dpi oc 5{p1)), 
ensuring that the transverse momentum integrated cross section is finite, 
while leaving problems at small finite transverse momenta. 

In the CCS approach |1] the resummed cross section for W production 
can be written as 

d(j 

— = ^o7(^^)^ 

t,W 

qq 

Xl ■ 



— = a^l{2^fY.\Vqq,\' / d%Upth)W^r{QAxi.X2 

dycL ptyv , J 



ZTT ZTX 



X2 



q2{x2,^^{b)) + — — — C^cq2{x2, Kb)) + — — — [Cg^g ® g2){x2, i^{b)) 

"1 



exp 



2 



iSgib,Q)+Sgib,Q)) 



(6) 

where the exponents in the Sudakov-like form factors read 



Sgib,Q)=Sgib,Q) 



djj 



n2 



2 

min 



In ( — ) + Bq{as{^Ji)) 



■ (7) 



The coefficient A and B can be expanded in the series of a^: 

A, = 2C. + i^^(...) + ..., 

/TT 

Bg = -3C^ + ^(...) + .... (8) 

The CSS formalism ^ leaves open the question of small b. Different pre- 
scriptions have been proposed to treat this region (see e.g.[7j). The lower 
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limit of the integral in Eq.Q is usually taken /i^j„ = (^) , where Cb = 2 
exp(— 7^;) ^ 1.1229. This prescription leads to a kink for the Sudakov form 
factor if Cb/b = Q. In order to allow a smooth dependence of the cross 
section on the gauge-boson transverse momentum and to guarantee that the 
lower limit is really lower than the upper limit, one could make the following 
replacement /i^^^ = (^)^ ^ (^)^ [1 + C^/{b'^Q^)]~^ (see e.g. 0). In order 
to guarantee that the scale of parton distribution does not take unphysically 
small values we shall use in addition the following prescription: 

/^'(&)=/^L„ + /^o, (9) 

where /Xq ^ ^qcd is the starting value for the QCD evolution. For illus- 
tration, in the present paper we shall use easy to handle the leading order 
parton distributions from Ref . |2Uj . 

In the formalism of unintegrated parton distributions the nonzero trans- 
verse momenta of gauge bosons are obtained already in the leading order. 
The invariant cross section for inclusive gauge boson production reads then 
as 

da _ W ST^ IT/ 12 r d'^Ki d?K2 



In the equation above the delta function assures the conservation of trans- 
verse momenta in the qq' fusion subprocess. The momentum fractions are 
calculated as xi^2 = exp(±y), where in contrast to the coUinear case 
M\Y is replaced by the transverse mass rrit^w Iii the case of the gauge boson 
production, the scale /i^ is taken as or M|. 

We note that formally, as far as matrix elements are considered, formula 
()1()|1 is the /cf-factorization counterpart of the collinear formula 0- 

Introducing unintegrated parton distributions in the space conjugated to 
the transverse momenta [TT ^ 

/g(a;, K^/i^) = ^ y exp (^inbj /g(x,6,/i^) (fb , (11) 
and taking the following representation of the 6 function 

5^(ki + K2 - Pt) = j d% exp (ki + K2 - Pt)b , (12) 

"'^In the present paper we shall use the notation / instead of / as in R,efs. |i;^l I14j to 
avoid confusion with antiquark distributions needed here. 
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the formula ()10|) can be written in the equivalent way 

/g/l(xi,6,/i2)/-,/2(x2,6,/i^) + /g7l(xi,6,/i2)/g/2(x2,&,yU^) • (13) 

In the next section we show and discuss the results obtained with the formula 
(fT!?j) . In the formulae for boson production \ Vqqi\'^ is replaced by Sqqi^{Vg + 

As already mentioned in the introduction, it is our intention here to use 
uPDFs /^^^^^(x,6,/i2) which fulfill b-space CCFM equations HUCg. How- 
ever, the perturbative solutions /^'^^^{x, b, iJ?) do not include nonperturba- 
tive effects such as, for instance, intrinsic momentum distribution of partons 
in colliding hadrons. In order to include such effects we propose to mod- 
ify the perturbative solution f^^^^ {x,h, ^'^) and write the modified parton 

distributions fq{x,b, fi"^) in the simple factorized form 

hix, b, ^^') = /f ^™(x, 6, ^^') ■ F^^{h) . (14) 

In the present study we shall use two different functional forms for the form 
factor 

F^ib) = F^^ib) = exp or exp (-^^ (15) 

identical for all species of partons. In Eq. (fT3j) 60 (or ^e) is the only free param- 
eter. In the next section we try to adjust this parameter to the experimental 
data on transverse momentum distribution of W^. 

We wish to note an almost identical structure of the b-space resummation 
formula ® and our formula (fT^ . Assuming a factorizable form for the 
nonperturbative form factor W^,^: 

W^YiQ^ ^ ^1' ^2) = i^,^^(Q, b, xi) ■ Ff^iQ, b, X2) , (16) 

the formulae © and (fT^ are identical as far as the formal structure is con- 
sidered, if we define effective b-space unintegrated (anti) quark distributions 

/f^^^(a;i,6,Q2) = ^Ff ^(Q, 6, xi) [xigi(xi, /i(&)) + ■••] exp Qs, (6, Q) 

/f «(X2,6,Q^) = ^Fi)^^(Q,6,X2) [x,q[{x„ fi{b)) + ...]exp Qs,-, {b,Q)(^7) 

The index SGR above stands for "soft-gluon resummation" . The fct-dependent 
unintegrated distributions of (anti) quarks corresponding to the b-space re- 
summation can be then obtained through the Fourier-Bessel transform 

/f ^^(a;, Q') = [ dbb UKb)f^''\x, b, Q'). (18) 
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With the simple Ansatz (fT^ for F^^ the whole dependence resides 
exclusively in the Sudakov-like form factor. 

3 Results 

For further reference in Figd we present results of the collinear NLO cal- 
culation. In this calculation the argument of the running strong coupling 
constant = and the scale argument is //^ = M^. In panel (a) we 
show the low-pt region, whereas in panel (b) the whole measured region. The 
\ow-pt problems are clearly visible. At the Fermilab energy the annihilation 
contribution dominates over the Compton one in the entire measured re- 
gion. The 2 — i> 2 NLO contribution underestimates the experimental data, 
especially in the interval 5 GeV < pt < 50 GeV. 

In our CCFM calculation we have taken /x^ = as the factorization 
scale in the b-space unintegrated parton distributions. In FigElwe present a 
two-dimensional map {y,Pt) of the invariant cross section. Large transverse 
momenta of W bosons are generated in the leading order calculations, which 
is due to transverse momenta of fusing partons contained in unintegrated 
(anti) quark distributions. 

In Figiniand Fig|3]we show transverse momentum distribution (integrated 
over rapidities) of in proton-antiproton collissions at Fermilab at W = 1.8 
TeV for the CCFM parton distributions (left panel) and standard soft-gluon 
resummation (right panel). The three curves in FigOlshow results obtained 
according to formula with different values of the Gaussian form factor 
parameter bo: bo = 0.5 GeV~^ (dashed), bo = 1.0 GeV~^ (solid) and bo = 2.0 
GeV~^. Similarly, the curves in Fig 0] correspond to fog = 0.5 GeV~^ (dashed), 
be = 1.0 GeV~^ (solid) and b^. = 2.0 GeV~^ in the exponential form factor. 
The results are overimposed on the DO collaboration data ^Sj measured at 
Fermilab. In order to convert the measured cross section {eu channel) to the 
cross section for W production we have taken the corresponding branching 
ratio BR{ev) =0.1073 jS]. For the Gaussian form factor (see FigO)) somewhat 
better agreement with the data is obtained with the CCS approach, while 
for the exponential form factor (see FigH} it is the opposite. One should 
remember, however, that the situation may change if 2 ^ 2 NLO processes, 
neglected here, are included to the CCFM formalism. In the region of large 
momenta those contributions are expected to be similar in size to the collinear 
NLO contributions shown in Fig^ (this is roughly the missing strength at pt 
= 20 GeV). In the CCS approach one somewhat arbitrarily matches the soft- 
gluon resummation result with the fixed order collinear result at the certain 
(chosen ad hoc) value of pt- We hope that in the formalism of unintegrated 
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parton distributions in the next-to-leading order approximation this can be 
done less arbitrarily. The DO data for W production are not precise enough 
to determine the parameter Bq (or fog) of the nonperturbative form factor 
with high precision. The low-pt DO data for production of Z boson should 
be better in this respect. 

In Fig0we show transverse momentum distribution (rapidity integrated) 
of Z bosons in proton-antiproton collisions at Fermilab at W = 1.8 TeV. As 
for the W production the theoretical results are compared with the DO exper- 
imental data pp. We have chosen the somewhat older data set P instead of 
the most recent one jB] because of the binning with slightly better resolution 
at small pt- In order to convert the measured cross section (e+e~ channel) to 
the cross section for Z production we have taken the corresponding branching 
ratio BR{e~^e~) =0.033632 The DO data allow better extraction of the 
parameter Bq of the nonperturbative form factor Fnp. By comparison to the 
data we find 60 = 0.5-1.0 GeV~^ (Gaussian) as giving the best description 
of the \ovf-pt Z^-pro duct ion data. 

Comparing FigH and FigsOl and |3] one observes that the leading-order 
CCFM approach generates more W's with 10 GeV < pt < 20 GeV than the 
standard next-to-leading order coUinear approach. The leading-order CCFM 
approach leaves some room at > 5 GeV. This is especially visible for the 
Gaussian form factor. In this context it would be very interesting to calculate 
the cross sections with the help of CCFM uPDF's up to the next-to-leading 
order. 

In the future, the gauge bosons will be measured also in proton-proton 
collisions at the BNL RHIC and at the CERN LHC. As an example of our 
predictions in FiglHlwe present the transverse momentum distribution of 
for three different center-of-mass energies W = 200, 500, 14000 GeV. In this 
calculation -2 < y < 2 and the parameter Bq was fixed at 1 GeV^^. The 
larger energy, the smaller-x region of (anti) quark distributions is probed. 
In principle, the form factor in ()15|1 can be not only a function of impact 
parameter b, but also of partonic momentum fraction x. Having precise data 
at different energies could help then to pin down the x dependence of the 
nonperturbative form factor. Other partonic processes could also be helpful 
in this respect. 

The precise determination of parton distributions is one of the most im- 
portant goals of high-energy physics. It was suggested in Ref . [22j that studies 
of the difference in the production of and W~ in proton-proton collisions 
could put further constraints on the d — u asymmetry as obtained from other 
sources, like deep-inelastic inclusive and exclusive scattering or Drell-Yan 
processes. In this illustrative leading order calculations, however, one could 
not address the problem of transverse momentum dependence of the W~ /W~^ 
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ratio. In contrast, it is very simple in our leading-order CCFM approach. In 
Figl?! we present the ratio as a function of W-boson transverse momentum 
at a typical RHIC energy W = 500 GeV for selected W-boson rapidities y=0 
(left panel) and y=± 1 (right panel). The solid curve corresponds to the cal- 
culation with the d — u asymmetry, while the dashed curve shows the result 
with d — u asymmetry switched off. In general, the ratio R = a{W~)/a{W~^) 
is smaller if the d — u asymmetry is included. The relative effect is larger for 
y=0. Fig. 13 illustrates therefore the potential to study the d — u asymmetry 
in W-boson production at RHIC. There is only a very weak dependence on 
W-boson transverse momentum. In order to illustrate the dependence on 
W-boson rapidity in FiglSlwe display also the two-dimensional maps {y,pt) 
of the asymmetry defined as 

aiW^)-aiW-) 
^ a{W+)+a{W-) ^ ' 

for asymmetric (left panel) and symmetric (right panel) antiquark distri- 
butions. The asymmetry so-defined is an asymmetric function under the 
operation y —>■ —y for proton-antiproton collisions and a symmetric one for 
proton-proton collisions. The STAR collaboration intends to measure the 
W-boson production in a few year perspective. It is not clear to us whether 
the map shown in FiglHl can be studied experimentally in the near future at 
the Brookhaven RHIC. 



4 Conclusions 

The zeroth order coUinear calculation leads to the 5-like distribution in the 
W-boson transverse momentum. The 1-st order calculation overestimates the 
experimental data at < 5 GeV and underestimates the data at pt > 5 GeV. 
In contrast, already the leading order CCFM approach leads to a reasonable 
transverse momentum distributions of W-bosons. The small transverse mo- 
menta pt < 10 GeV are sensitive to the nonperturbative physics embodied in 
our calculation in the extra b-dependent form factors. The parameter(s) of 
the form factor can be adjusted to the experimental low-pf distributions pro- 
vided the data are sufficiently precise. The simple leading order calculation 
leaves some room for next-to-leading order contributions at pt > 5-10 GeV. 

In the present study we have used very simple functional forms of the 
nonperturbative form factors responsible for intrinsic momentum distribution 
of partons in colliding hadrons. It is not obvious a priori if these forms 
are sufficient to describe the variety of data including Drell-Yan dilepton 
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production, prompt photon production, heavy quark pair production, etc. 
This will be a subject of the forthcoming research. 

We have studied the dependence of the cross section ratio W~ /W^ in 
proton-proton collisions on transverse momentum and rapidity of the W- 
boson. We find a potential to study the d — u asymmetry at BNL RHIC 
quantitatively. 

There are many other possible ways to study the unintegrated parton 
distributions further. For instance at RHIC energy the topics to study can 
be jet or/and particle production 

Acknowledgments This paper was completed after Jan Kwiecihski passed 
away. One of us (A.S.) is indebted to Krzysztof Golec-Biernat for a discus- 
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Figure 1: Transverse momentum distribution ofW^ + at W = 1.8 TeV. 
The curves represent the full first-order result (solid) and the annihilation 
(dashed) and Compton (dotted) components separately. The experimental 
data are taken from Ref-l^. 
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Figure 2: Two-dimensional map of the cross section for W~^ + W~ production 
in proton- antiproton collisions at W — 1.8 TeV. In this calculation bg = 1 
GeV~^ ( Gaussian form factor). 
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Figure 3: Transverse momentum distribution ofW^ + W~ at W = 1.8 TeV 
within the formalism of the CCFM uPDF's (left panel) and soft-gluon re- 
summation (right panel) with Caussian nonperturbative form factor. The 
three curves correspond to different values of the form factor parameter bo as 
explained in the text. The experimental data are taken from Ref.^. 
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Figure 4: Transverse momentum distribution ofW~^ + iy~ at W = 1.8 TeV 
within the formalism of the CCFM uPDF's (left panel) and soft-gluon re- 
summation (right panel) with exponential nonperturbative form factor. The 
three curves correspond to different values of the form factor parameter he as 
explained in the text. The experimental data are taken from Ref.'^. 
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Figure 5: Transverse momentum distribution of at W = 1.8 TeV. The 
three curves correspond to different values of the Gaussian form factor pa- 
rameter Bq as explained in the text. The experimental data are taken from 

Ref.m- 
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Figure 6: Transverse momentum distribution of W in proton-proton col- 
lisions for three different energies. In this calculation -2 < y < 2 and the 
parameter of the Gaussian form factor hQ — 1 GeV~^. 
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Figure 7: T/ie raiio cr(Vr~)/(T(Vr"'") /or proton-proton collisions at RHIC 
energy W—500 GeV for y—0 (left panel) and y—±l (right panel). In this 

calculation the parameter bo of the nonperturbative Gaussian form factor was 
fixed for b^ = 1 GeV^^ . The solid line corresponds to calculation with d — u 
asymmetry while the dashed line to calculation with symmetric sea. 
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Figure 8: Two-dimensional maps of the asymmetry A± for proton-proton 

collisions at RHIC energies W=500 GeV for two cases: in panel (a) with 
d — u asymmetry and in panel (b) without the asymmetry. In this calculation 
the parameter bo of the nonperturbative Gaussian form factor was fixed at bo 
- 1 GeV-\ 
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